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Abstract. The Hall conductivity in the mixed state of a clean (i > ^o) type-II s-wave 
superconductor is determined from a microscopic calculation within a quasiclassical 
approximation. We find that below the superconducting transition the contribution 
to the transverse conductivity due to dynamical fluctuations of the order parameter 
is compensated by the modification of the quasiparticle contribution. In this regime 
the nonlinear behaviour of the Hall angle is governed by the change in the effective 
quasiparticle scattering rate due to the reduction in the density of states at the Fermi 
level. The connection with experimental results is discussed. 
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The Hall effect in the mixed state of type-II superconductors has remained a 
theoretical puzzle for almost thirty years Q. The existing phenomenological || [J 
theories predict that the Hall angle in the flux-flow regime is either identical to that in 
the normal state |2j or constant ||. Theories which make use of the time-dependent 
Ginzburg-Landau equations (TDGL) also find no modification of the Hall conductivity 
in the superconducting state [|IJ. These predictions are at variance with the strongly 
nonlinear behaviour (as a function of magnetic field) found in experiments performed 
on both low-T c materials P, ^] and high-T c cuprates EJ. For dirty superconductors 
(I <C £o> where I is the mean free path and £o is the superconducting coherence length), 
transport coefficients can be determined from microscopic theory by a straightforward 
expansion in powers of the order parameter, A. The results of such a calculation for 
the transverse resistivity @ explain qualitatively the sharp increase in the Hall angle 
observed in experiment (although, to our knowledge, no systematic comparison has 
been made), and provide the physical basis for a generalized TDGL approach, in which 
the relaxation rate is assumed to be complex, rather than purely real, to allow for a 
modification of the transverse transport coefficients |], |§- The small parameter in the 
expansion of the microscopic equations is proportional to both the order parameter and 
the mean free path, therefore, it is not small in the clean (/ 3> £o) limit. In this regime a 
straightforward expansion is not possible and the TDGL equations are not applicable ||, 
so that an alternative approach is needed for the calculation of the transverse transport 
coefficients. 

In this Letter we give the results of a calculation of the Hall conductivity of a 
clean s-state superconductor in the mixed state near the upper critical field, H c2 , which 
uses a quasiclassical approximation to the microscopic theory. We made this choice as 
both the normal state and superconducting properties of the low-T c compounds are well 
known, and comparison between theory and experiment is fraught with less ambiguity. 
The quasiclassical approach []T0j| has been applied successfully in the past to study 
transport phenomena in superfluids |l l]] and superconductors [12] and more recently to 
the unconventional superconductors fl3| . The central quantity in this method is the 
single-particle matrix Green's function G integrated over the quasiparticle energy 

g(s, R; u n , u n <) = ^| = J ^G{p, R; u n , u n >); (1) 

here the ui n = 2irT(n + |) are fermionic Matsubara frequencies, s is the normalized 
parameterization of the Fermi surface, R is the center of mass coordinate and p is 
the relative momentum. Since the Green's function is strongly peaked at the Fermi 
momentum pf, which is normally far larger than any other momenta in the problem, 
slower varying quantities such as the self-energy and external potential can be expanded 
around their values at the Fermi surface. The result of such an expansion [O] in the 
small parameter l/^kf ~ A/e/ is a set of transport-like equations for the quasiclassical 
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propagator g. We have generalized these equations to include terms responsible for the 
Hall Effect in a charged superfluid. Technical details of the derivation will be reported 
elsewhere |b|], here we use the equations to determine the transverse dc-conductivity 



We use linear response theory in the vector potential A(co>o) describing a constant 
electric field E = E5t. The magnetic field H = Hz, chosen parallel to the z-axis, 
is described by the vector potential *4(R) = Hxy. We consider a spherical Fermi 
surface and use the Born approximation for s-wave impurity scattering characterized by 
a collision rate r -1 . The spatial dependence of the order parameter is modeled by the 
periodic Abrikosov vortex lattice 

A(R) = Y.C ky e^ y %{x - A%). (2) 

by 

Here $o( a; ) is the lowest energy eigenfunction of the linearized Ginzburg-Landau 
equation (i.e. the eigenfunction of a harmonic oscillator with Cooper pair mass M = 2m 
and frequency u c ) and A 2 = (2eH)^ 1 is the magnetic length ( A ~ £ for fields H ~ H c2 ). 
This approach is appropriate provided that the broadening of the levels in the vortex 
core is large compared to their spacing 1/r ^> A 2 /e/, it breaks down in the superclean 
regime (cf. jnj). In the clean limit the finite lifetime is accounted for by replacing 
uj n by ui n = u) n + (2r) _1 (g(u; n )) (angular brackets denote an average over the Fermi 
surface); corrections to the order parameter due to impurity renormalization are of the 
order 0(A/l) and can be ignored. The equations for the unperturbed functions / and 
g and the linear, in A, corrections to the propagator f\ and g\ are p 



[2u n + V/ (V - 2ieA)}f = 2iAg (3) 
[2Q n + v/(V - 2ieA)]fx = zev/A(/ + /(-)) + iA( 9l - g x ) + iA^ + <?(-)) (4) 
(iu + iu) c ^-){ 9l - 9l ) = 2ev f A(g - g(-)) + A*(f - /(-)) + A x (/t _ /t(_)) ( 5 ) 



+ (2r)- X ((/!)(/ - /(-)) + (h)(P - /+(-))) 



+2 dA*dfl , n dAdf x 



<9R dp\\ dHdp\\ 

where Ai is the change in the order parameter induced by the electric field. There 
are corresponding equations for f< and g. In these equations a z is the Pauli matrix, 
the Fermi velocity vj(s) = v (sin 9 cos0, sin 9 sin 0, cos#), p|| is the component of the 
momentum parallel to the Fermi surface, uj c = eH/mc is the cyclotron frequency, and the 
angular brackets denote an average over the Fermi surface. We have used the shorthand 
notations uj n — = uj n —ujQ andg = g(& n ), g{—) = g(u n —), and introduced 2Vt n = uj n +uj n — 
, uj = oj n — ijj n —. The superconducting order parameter satisfies the usual self- 



4 



consistency condition A(R) = gN(0)nT J2 n I d 2 sf(s, R; u n , u n ), where g is the coupling 
constant and iV(0) is the density of states at the Fermi surface. The unperturbed 
function g obeys the normalization conditions g + g = and g 2 — ff> = — 1 ||10|| . We 
chose to write an equation for g\ — g~\ since the transport current is given by |T(| 

j = ivreiV(0)T^ f d 2 s y f (s)(g 1 (s,u n ) -g~i(s,u n )). (6) 

and the functions fx, f\ depend on gx, gi in this combination only In equations (3)- (5) 
we have omitted terms whose contribution to the conductivity is of order 0(A/l) smaller 
than that of leading order terms. Equation (3) is the well-known static Eilenberger 
equation [|l(J] and equations (3)- (5) contain all the terms relevant to the Hall effect for a 
clean type-II superconductor in the high-field region. It is clear from equation (5) that 
there are several distinct contributions to g± — g\ (and therefore the current). The first, 
the quasiparticle contribution, depends on the unperturbed function g, while the second, 
proportional to A 1; is due to the dynamical fluctuations of the order parameter induced 
by the perturbing electric field. The third term in this equation describes the additional 
scattering of quasiparticles off these dynamic fluctuations; it has the same origin as 
the Thompson diagram in the analysis of transport in dirty superconductors [pjj. The 
fourth term describes how as vortices move and are deformed by the transport current, 
the resulting gradients of the order parameter act as driving forces in the transport- 
like equations. Finally, as the renormalization of the frequency u Q in equation(5) 
depends on the angular average of the quasiparticle propagator, which changes in the 
superconducting state, the effective transport mean free path is modified. 

To proceed we approximate the diagonal part of the quasiclassical propagator by 



its spatial average [I8|. Since the electromagnetic fields in a superconductor vary over 
distances of the order of the penetration depth A, this is a very good approximation 
in the London limit k = A/£o ^ 1; even for compounds with moderate values of k it 
remains valid for a wide field range below H c2 . In all of the following g will stand for 
the averaged distribution function. We now solve equation(3) with the normalization 
condition to determine the unperturbed functions / and g. We then determine f\ 
and Ai to leading order in (AA/t>) 2 C 1 by solving equation (4) together with the 
self-consistency condition for the order parameter. To accomplish this program the 
expression (2u n + Vf ■ (V — 2ieA))~ l A has to be evaluated. To do this we exploit 
the oscillatory character of the Abrikosov solution and introduce raising and lowering 
operators a = (A/y/2)\y x + i(V y - 2ieHx)} and a) = -(A/V2)[V X - i(V y - 2ieHx)} 
obeying the usual bosonic commutation relations [a,a^] = 1. If the ground state 
equation (0) is denoted by |0), the higher eigenstates (modes of the order parameter) 
are generated by the standard procedure cr |n) = \/n + l\n + 1). Then to make use 
of the properties of these operators the operator vj(V — 2ieA) can be rewritten as 

and the result of its acting on any state \n) evaluated 
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explicitly |b|]. Using this approach we are able to determine the unperturbed functions 



g = -i&gn(u r . 



2AA 2 . 2iuj n Asga(uj n ) , 



-1/2 



, .. . . (7) 

" v sin v sin i 

where A is the spatial average of the order parameter and W(z) = e~ z2 erk(—iz), and 

f o- V^ A ^ 1 (- i \ m ( ( \\m+l irrubuHm) / 2i£ n Asgn(^ n ) \ 

The expression for g reproduces correctly the gapped BCS-like function for quasiparticles 
traveling parallel to the magnetic field, while describing gapless behaviour in all other 
directions. A similar expression for g has been obtained by Pesch fT9"| . Since / is a 



Fourier series in <p, the mode with m = will couple to a scalar potential, the mode 
with m = 1 to a transverse potential etc. Then we find A 1 = (ieAA\^2)[(l — iur)^ 1 + 
uj c tujt\\1) (here to is the real external frequency). Using this value to determine the 
correction to anomalous propagator j\ and calculate the current from equation (Q) we 
obtain the longitudinal and transverse conductivities up to order (AA/v) 2 . 

If we use the usual notation a n for the normal state conductivity a n = N(0)e 2 v 2 r/3, 
the enhancement of the transverse current due to Lorentz force driven fluctuations of 
the order parameter 

ag = Qa n uj c r(AA/v) 2 (9) 

is exactly compensated by the modification of the quasiparticle Hall current due to 
additional scattering off the vortex lattice 

Aaf y = -Qa n u c T{AA/v) 2 . (10) 

Similarly, the positive contribution to the transverse conductivity due to the forces 
generated by gradient of the excited mode of the order parameter 

a 9 J y = 3a n (A 2 r/E f ) = 12a n Lu c r(AA/v) 2 (11) 

is cancelled by the additional scattering introduced by the deformed and moving vortex 
lattice 

a™ = -\2a n u c r{AA/v) 2 . (12) 

As a result, the behaviour of the transverse conductivity a xy is determined solely by the 
effect of the modification of the effective elastic scattering time T e ff on the leading order 
quasiparticle contribution. For the dc conductivity this change is due to the decrease in 
the number of states at the Fermi surface available for scattering as the superconducting 
gap opens. We find, in agreement with the result of Pesch [IB3] for the density of states, 
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that the increase in the relaxation time is a non-analytic function of the small parameter 
(AA/u) 



'eff 



1+4 



/AA\ 2 /AA 



/AA 



and, up to order (AA/v) 2 the transverse conductivity is given by 



a, 



xy 



~N(0)e 2 v 2 T eff (u; c T eff ) 



a n u c r 



1 + 4 



/AA 



I" 



log 



2v 2 
AW 



The longitudinal conductivity is obtained in a similar way 



Or, 



- i 



(13) 



(14) 



(15) 



In the high-field region the square of the order parameter is linear in the applied magnetic 
field and is given by pUfl 



A 2 



H c2 — H 



TcN(0)f3 A (2K 2 2 -l) 



H, 



c2 



T dH c2 
2 dT 



(16) 



In figure 1 we show the qualitative behaviour of function of magnetic 

field, plotted using parameter values for Nb. The transverse conductivity is enhanced 
below the upper critical field and has negative curvature in the high field region. While 
the transverse conductivity is proportional to the square of the scattering time, the 
Hall angle tan## = u xy /u xx is only linearly dependent on the scattering time and 
the corresponding nonlinear dependence on magnetic field is weaker, as can be seen 
in figure 2. Finally, as the transverse resistivity p xy ~ &xy/&1 x is independent of the 
effective scattering time, upon entering the superconducting state it remains linear in 
magnetic field with the same slope as in the normal metal. This behaviour is to be 
contrasted with that of Bardeen-Stephen model H, where the resistivity is modified but 
the Hall angle obeys the same linear law as in the normal state. On the other hand, the 
Nozieres-Vinen theory @, which predicts that the Hall angle should be constant in the 
flux-flow regime below H c2 at variance with the result of this work, also finds that the 
transverse resistivity is identical to that of the normal state, although the individual 
components of the conductivity tensor are quite different from those found here. 

A comparison can be made with the experimental data of Fiory and Serin || 
on high purity Nb. These experiments find a transverse resistivity in the flux-flow 
regime which is linear in the applied magnetic field over a wide range of fields below 
H c2 . The Hall angle, however, flattens or even increases above its value at H c2 before 
decreasing at lower fields. These results are more suggestive of the behaviour given here 
than the original interpretation given in terms of the Nozieres-Vinen theory. Also, the 
longitudinal resistivity found in |J has a distinct increase in slope just below the upper 
critical field, which is consistent with the logarithmic behaviour given by equation flT5|). 
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Such comparisons are, of course, only qualitative, and more experimental work is needed 
to make a detailed comparison with the theory. 

To conclude, we give the results of a microscopic calculation of the Hall resistivity 
of a clean type-II s-state superconductor in the high- field limit. We find that the field 
dependence of the Hall conductivity in the high field regime, which is non-analytic, 
is entirely due to the change in the density of quasiparticle states at the Fermi level 
in the superconducting state. At the same time we find that the field dependence 
of the transverse resistivity below the upper critical field remains unchanged. These 
results are in qualitative agreement with the experimentally observed behaviour. The 
approach developed here can be generalized to superconductors with unconventional 
order parameter symmetry, this work is now in progress. 

One of us (A H) would like to thank A T Dorsey, D Rainer, and K Scharnberg for 
important discussions. This research was supported in part by the National Science 
Foundation through Grant No DMR9008239. 

References 

[1] For a recent review see Dorsey A T 1992 Phys. Rev. B 46 8376 

[2] Bardeen J and Stephen M J 1965 Phys. Rev. 140 A1197 

[3] Nozieres P and Vincn W F 1966 Philos. Mag. 14, 667 

[4] Hagcn S J et al 1990 Phys. Rev. B 41 11630 

[5] Fiory A T and Serin B 1968 Phys. Rev. Lett. 21 359 

[6] Harris J M et al 1996 J. Low Temp. Phys. 105 877 

[7] Ebisawa H 1972 J. Low Temp. Phys. 9 11 

[8] Maki K 1969 Phys. Rev. Lett. 23 1223 

[9] Houghton A and Maki K 1971 Phys. Rev. B 3 1625 

[10] Eilenbergcr G Z. Phys. 1968 214 195 

Larkin A I and Ovchinnikov Yu N 1968 Zh. Exp. Teor. Fiz. 55 2262 [Sov. Phys. JETP 1969 28 
1200] 

[11] Serene J W and Rainer D Phys. Rep. 1983 101 221 and references therein. 

[12] Ovchinnikov Yu N Zh. Eks. Teor. Fiz. 1974 66 1100 [Sov. Phys. JETP 1974 39 538] 

[13] Graf M J, Palumbo M, Rainer D and Sauls J A, 1995 Phys. Rev. B 52 10588 

[14] Houghton A and Vekhter I (to be published). 

[15] Kopnin N B 1994 J. Low Temp. Phys. 97 157 

[16] Larkin A I and Ovchinnikov Yu N 1995 Phys. Rev. B 51 5965 

[17] Thompson R 1970 Phys. Rev. B 1970 1 327 

[18] Brandt U, Pesch W and Tcwordt L 1967 Z. Phys. 201 209 

[19] Pesch W 1975 Z. Phys. B 21 263 

[20] Maki K and Tsuzuki T 1965 Phys. Rev. 139 A868 



8 

Figure captions 



Figure 1. Hall conductivity as a function of the reduced magnetic field 



Figure 2. Hall angle as a function of the reduced magnetic field 



